Atmospheric 
Chemistry 
and Physics 


Atmos. Chem. Phys., 11, 6245-6263, 2011 
www.atmos-chem-phys. net/1 1/6245/201 1/ 
doi:10.5194/acp-l 1-6245-2011 
© Author(s) 2011. CC Attribution 3.0 License. 



Simultaneous retrieval of aerosol and cloud properties during the 
MILAGRO field campaign 

K. Knobelspiesse 1 , B. Cairns 2 , J. Redemann 3 ’ 4 , R. W. Bergstrom 3 , and A. Stohl 5 

1 NASA Postdoctoral Program Fellow, NASA Goddard Institute for Space Studies, New York, New York, USA 

2 NASA Goddard Institute for Space Studies, New York, New York, USA 

3 Bay Area Environmental Research Institute, Sonoma, California, USA 

4 NASA Ames Research Center, Moffett Field, California, USA 

5 Norwegian Institute for Air Research, Kjeller, Norway 

Received: 28 January 201 1 - Published in Atmos. Chem. Phys. Discuss.: 23 February 201 1 
Revised: 3 June 201 1 - Accepted: 14 June 201 1 - Published: 1 July 201 1 


Abstract. Estimation of Direct Climate Forcing (DCF) due 
to aerosols in cloudy areas has historically been a diffi- 
cult task, mainly because of a lack of appropriate measure- 
ments. Recently, passive remote sensing instruments have 
been developed that have the potential to retrieve both cloud 
and aerosol properties using polarimetric, multiple view an- 
gle, and multi spectral observations, and therefore determine 
DCF from aerosols above clouds. One such instrument is 
the Research Scanning Polarimeter (RSP), an airborne pro- 
totype of a sensor on the NASA Glory satellite, which un- 
fortunately failed to reach orbit during its launch in March 
of 2011. In the spring of 2006, the RSP was deployed on 
an aircraft based in Veracruz, Mexico, as part of the Megac- 
ity Initiative: Local and Global Research Observations (MI- 
LAGRO) field campaign. On 13 March, the RSP over flew 
an aerosol layer lofted above a low altitude marine stratocu- 
mulus cloud close to shore in the Gulf of Mexico. We in- 
vestigate the feasibility of retrieving aerosol properties over 
clouds using these data. Our approach is to first determine 
cloud droplet size distribution using the angular location of 
the cloud bow and other features in the polarized reflectance. 
The selected cloud was then used in a multiple scattering ra- 
diative transfer model optimization to determine the aerosol 
optical properties and fine tune the cloud size distribution. In 
this scene, we were able to retrieve aerosol optical depth, the 
fine mode aerosol size distribution parameters and the cloud 
droplet size distribution parameters to a degree of accuracy 
required for climate modeling. This required assumptions 
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about the aerosol vertical distribution and the optical proper- 
ties of the coarse aerosol size mode. A sensitivity study was 
also performed to place this study in the context of future sys- 
tematic scanning polarimeter observations, which found that 
the aerosol complex refractive index can also be observed ac- 
curately if the aerosol optical depth is larger than roughly 0.8 
at a wavelength of (0.555 pm). 


1 Introduction 

The radiative effects of atmospheric aerosols (suspended par- 
ticles) are among the least certain components of global cli- 
mate models (IPCC, 2007). Even the most simple aerosol 
radiative effect. Direct Climate Forcing (DCF), is difficult to 
model, mainly due to the heterogeneity of aerosol sources, 
evolution, sinks, and radiative properties, and the difficulties 
of global observation (Mishchenko et al., 2009). DCF is de- 

fined as the change in net irradiance at the top of the atmo- 

sphere from the scattering and absorption of anthropogenic 

aerosols, and it neglects more complicated impacts due to 
aerosol induced changes in cloud properties or the atmo- 
spheric temperature vertical profile. Current modeling ca- 
pability can be found in the results of the AeroCom model 
intercomparison effort. For example, Schulz et al. (2006) 
found a DCF range of +0.04 to —0.41 (W m -2 ). A sensi- 
tivity analysis in this work found that DCF is strongly af- 
fected by aerosol forcing efficiency (the radiative forcing of a 
unit optical depth). Unfortunately, many factors that change 
aerosol forcing efficiency, such as particle size, absorption, 
and refractive index, are rarely observed to the level of accu- 
racy required by climate models (Mishchenko et al., 2007b). 
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Furthermore, most aerosol remote sensing algorithms with 
passive instruments are successful only in cloudless condi- 
tions. Efforts to observe aerosols in the presence of clouds 
are still in their infancy, and are often intertwined with 
aerosol optical property assumptions. 

The object of this study is the retrieval of aerosol prop- 
erties when they are suspended above liquid boundary layer 
clouds. Aerosols above clouds (henceforth abbreviated as 
AAC) are a potentially important component of positive 
DCF, since absorbing AAC can significantly reduce a high 
underlying cloud albedo and therefore alter forcing (Hay- 
wood et al., 1997). However, reliable estimates of the global 
prevalence and anthropogenic component of this type of 
forcing are limited by a lack of appropriate observations, es- 
pecially since DCF is strongly dependent on aerosol optical 
properties. Furthermore, AAC can interfere with the ability 
of passive instruments to accurately determine cloud opti- 
cal properties (Coddington et al., 2010). Several approaches 
have been developed recently to observe AAC, but they are 
limited in their ability to distinguish aerosol types because 
of the significant assumptions required by their retrieval al- 
gorithms. Chand et al. (2008) used the active observations 
of the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite 
Observations (CALIPSO) instrument to determine AAC op- 
tical depth at two wavelengths. The ratio of the optical depth 
spectral pair suggests the aerosol particle size. This method 
is therefore somewhat limited in its ability to determine DCF 
from AAC, although these type of data show promise in 
their potential for combined retrievals with passive remote 
sensing data. Another method uses passive spectrometer ob- 
servations at ultra-violet wavelengths from instruments such 
as the Scanning Imaging Absorption Spectrometer for At- 
mospheric Chartography (SCIAMACHY). De Graaf et al. 
(2007) fit simulations of biomass burning (smoke) AAC to 
observations from the SCIAMACHY instrument to deter- 
mine the aerosol total and absorbing optical thickness. An al- 
ternate technique, which we apply here, is to use multi-angle, 
multi-spectral, passive observations of polarized reflectance 
to simultaneously determine cloud and AAC optical prop- 
erties. Waquet et al. (2009b) demonstrated this technique 
using data from the Moderate Resolution Imaging Spectro- 
radiometer (MODIS) and the Polarization and Anisotropy 
of Reflectances for Atmospheric Sciences Coupled with Ob- 
servations from a Lidar (PARASOL) instruments. Both in- 
struments were a part of the so called “A-train” polar orbit 
until recently (January 2010), when the PARASOL instru- 
ment was moved to another orbit with the consumption of the 
fuel it needed to safely maintain its position. Waquet et al. 
(2009b) combined MODIS retrievals of cloud top height with 
PARASOL polarized observations at a variety of scattering 
angles to determine the aerosol optical thickness of biomass 
burning AAC in the South Atlantic Ocean. This required as- 
sumptions about the aerosol size and refractive index, since 
a single PARASOL band (at 0.865(pm)) and a single scat- 
tering model were used. PARASOL also has a somewhat 


coarse angular resolution, which means observations over a 
large area must be collected in order to observe the cloud 
bow (which is similar to a rainbow, but is due to scattering 
from cloud, and not rain, droplets) with enough angular sam- 
pling to accurately constrain cloud optical properties. Nev- 
ertheless, we believe this approach has potential for future 
instruments that have a higher angular resolution and more 
available spectral bands than PARASOL. 

The Aerosol Polarimetery Sensor (APS) is a scanning po- 
larimeter was to be an instrument on the NASA Glory satel- 
lite (Mishchenko et al., 2007a). Unfortunately, Glory failed 
to reach orbit during its launch in March of 2011. Despite 
this setback, deployment of the APS airborne prototype, the 
Research Scanning Polarimeter (RSP) continues, and so does 
the analysis of the capabilities of this class of instrument. 
Both instruments are is intended to retrieve aerosol and cloud 
optical properties to a high degree of accuracy. They do so 
by gathering a large amount of information in each scene: so- 
lar reflected linear polarization at about 250 viewing angles 
in nine channels at visible and near infra-red wavelengths. 
Cloud and aerosol properties are retrieved by matching a 
doubling and adding radiative transfer model (Hansen and 
Travis, 1974) to observations. The Research Scanning Po- 
larimeter (RSP), was developed to test APS observational ca- 
pabilities (Cairns et al., 1997; Cairns, 2003). APS and RSP 
have very similar characteristics, where the latter has fewer 
view zenith angles (152) and slightly different center wave- 
lengths for some bands. The RSP has flown on a variety of 
aircraft in many field campaigns. Field campaign data have 
been used to validate the ability of RSP/ APS type instru- 
ments to observe aerosols over the ocean (Chowdhary et al., 
2001; Chowdhary et al., 2002; Chowdhary et al., 2005a), 
aerosols over land (Elias et al., 2004; Waquet et al., 2009a), 
cloud optical properties (Chowdhary et al., 2005b), surface 
reflectance properties (Knobelspiesse et al., 2008; Litvinov 
et al., 2010), and of extremely optically thick plumes where 
the surface is obscured (Knobelspiesse et al., 201 1). 

In the spring of 2006, the RSP participated in a large 
field campaign called Megacity Initiative: Local and Global 
Research Observations (MILAGRO). The goal of MILA- 
GRO was to examine the behavior, transport, and evolu- 
tion of emissions from the Mexico City Metropolitan Area 
(MCMA). More than 450 scientists from 150 institutions 
participated, involving instrumentation onboard seven air- 
craft, multiple ground sites, and the coordinated use of at- 
mospheric models and orbital observations (Molina et al., 
2010). NASA’s participation in MILAGRO was called the 
Intercontinental Chemical Transport Experiment-B (INTEX- 
B), for which the RSP was deployed on the Sky Research, 
Inc., Jetstream-31 (J-31) aircraft. The J-31 was based in 
Veracruz, Mexico for three weeks in February and March, 
and performed a total of thirteen successful research flights. 
In addition to the RSP, several other instruments were part 
of the J-31 payload. This included the Ames Airborne 
Tracking Sunphotometer (AATS-14) (Redemann et al., 2009; 
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Livingston et al., 2009), the Solar Spectral Flux Radiometer 
(SSFR) (Pilewskie et al., 2003; Bergstrom et al., 2010), the 
Cloud Absorption Radiometer (CAR) (Gatebe et al., 2010) 
and a Position and Orientation System (POS) and various 
meteorological sensors. 

Aerosols encountered during MILAGRO were a compli- 
cated mix including urban/industrial particles from Mexico 
City (Doran et al., 2007; Lewandowski et al., 2010; Marley 
et al., 2009; Paredes-Miranda et al., 2009; Rogers et al., 2009 
and others referenced in Molina et al. (2010)), smoke from 
local agricultural fires, and even volcanic sulfates (Grutter 
et al., 2008, de Foy et al., 2009). In the first two weeks of 
March, 2006, transport from the MCMA was weak but gen- 
erally to the east, allowing aerosols to mix and age before 
they were transported over the Gulf of Mexico. On the morn- 
ing of March 13th, a layer of low altitude (probably marine 
stratocumulus) clouds formed over the Gulf coast of Mex- 
ico between the cities of Veracruz and Tamiahua (approxi- 
mately 96.75° W and 20.25° N). The cloud slowly dissipated 
throughout the day, but was overflown by the J-31 aircraft, 
which later descended to the surface in a recently cloud free 
region. Observations by the SSFR and AATS-14 instruments 
on the J-31, along with aerosol transport model simulations, 
suggest that a layer of mixed aerosols were lofted above the 
cloud. This is therefore an ideal test scene for the ability of 
sensors such as RSP and APS to retrieve optical properties of 
AAC. 

In the next section of this paper (2), we describe our 
methodology for the retrieval of AAC optical properties. 
Next (Sect. 3), we perform a short sensitivity study to help 
us construct an appropriate optimization strategy for AAC 
scenes. We then examine the retrieval results, and compare 
them with observations from other instruments on that day 
(Sect. 4). We also investigate the optimization sensitivity to 
forward model assumptions. We then discuss the implica- 
tions of our retrieval and sensitivity studies on the ability of 
RSP and APS to resolve AAC optical properties in the con- 
text of climate model requirements (Sect. 5), followed by a 
brief conclusion. 


2 Method 

One benefit of multi-angle polarized measurements is that 
the various atmospheric components in a scene modify dif- 
ferent angular portions of the observation. For our scene, 
where aerosols overlay a warm phase cloud, the cloud wa- 
ter droplet size distribution determines the angular location 
of the cloud bow (which typically exists at scattering angles 
somewhere around 142°) and the magnitude of secondary 
bows (Breon and Goloub, 1998). Aerosol properties such as 
optical depth, size distribution, and refractive index dominate 
the signal at side scattering angles less than 120° (Goloub 
et al., 2000). Rayleigh scattering also contributes at these 
angles, which means multi-angle polarization observations 


555nm 




Fig. 1. The sensitivity of multi-angle polarimetry to cloud and 
aerosol optical properties is demonstrated in this figure. Black lines 
are ( Rq ) in the scattering plane for a simulated cloud (type A) 
with an optical thickness of 20.0 and a vertically uniform droplet 
size distribution of ( r e =6.25 pm), (v e = 0.75) from the ground to 
480 m. The top panel is reflectance at 0.555 pm, while the bot- 
tom panel is reflectance at 0.865 pm. Blue lines are reflectance 
from the cloud when the top has been raised by 500 m. Green 
lines are the reflectance of a cloud containing different droplet sizes, 
({r e = 10 pm), (v e — 0.05)). Red and magenta lines indicate the re- 
flectance of a cloud with aerosols above, the former for “Mexico 
City” type urban aerosols with (r (0.555) = 0.2) and the latter for 
“African Savann” biomass burning aerosols with an optical thick- 
ness of (r(0.555) = 0.4) from Dubovik et al. (2002). All scenes are 
simulated with a solar zenith angle of (Q s =45°) and a relative az- 
imuth angle of (cf> = 45°). 

can be used to determine cloud top heights in the absence of 
aerosols (Goloub et al., 1994). When aerosols are present, 
retrieval of their optical properties is helped by external in- 
formation about the cloud top height. Figure 1 illustrates the 
available information in a synthetic AAC scene. This is ex- 
pressed by (Rq), which is defined below in Eq. (1). Note 
how the angular location of the cloud bow at high scattering 
angles changes with cloud droplet size distribution. Cloud 
top height has an impact at smaller scattering angles and the 
shorter wavelength, while AAC modify the magnitude of the 
cloud bow peak and alter the polarized reflectance at scatter- 
ing angles less than (120°) for all wavelengths. 
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In this paper, we test the capability of RSP to retrieve 
AAC optical properties. To do this, we start by identify- 
ing the cloud and aerosol vertical distribution using obser- 
vations from other instruments on the J-31 aircraft. These 
are used to determine layer heights in a plane parallel radia- 
tive transfer model, which we iteratively modify to match to 
our observations using the Levenberg-Marquardt technique. 
The aerosol and cloud optical properties that create the best 
match are the “retrieved” values. Initial optical properties 
must be close to the retrieved solution and therefore selected 
with care. We find cloud initial values by matching obser- 
vations of the cloud bow angular location and width to sin- 
gle scattering simulations (scattering phase functions com- 
puted with Mie theory). Aerosol initial values were selected 
from the “Urban-industrial and mixed; Mexico City” class 
described in Dubovik et al. (2002). This climatology was 
derived from a year of ground based measurements by the 
Aerosol Robotic Network (AERONET). The aerosol optical 
depth initial value comes from observations by the AATS-14 
instrument during a spiral to the surface in a cloud free region 
near our scene. 

The RSP observes the first three components of the Stokes 
polarization vector ((/), (Q) and (£/)), meaning it is sen- 
sitive to both linearly polarized and total radiance (Hansen 
and Travis, 1974). For remote sensing, it is often useful to 
express the polarization components of the Stokes vector in 
terms of reflectance. This takes the exo-atmospheric radiance 
into account, and is calculated: 


Inn 


Rl — F„cos6, 

R Q = 

Ru = 


Qx’o 
F„c osftj 

UnrZ 


F„cos0 s 


( 1 ) 


where ( F a ) is the annual average exo-atmospheric irradiance 
((W/m 2 )), ( r a ) is the solar distance in AU (thus compensat- 
ing for solar distance deviation from average throughout the 
year), and (9 S ) is the solar zenith angle. The RSP instrument 
has a high (0.2 %) absolute accuracy for (Rq) and (Ru) rela- 
tive to (Ri), as the same detectors are used to measure inten- 
sity and linear polarization, and because calibration of their 
relative gain is performed with each scan. (Ru) typically be- 
comes very small when ( Q) and (U) are defined with respect 
to the scattering plane (containing the solar illumination and 
observation vectors), so we subsequently use (Rq) defined 
in this plane as a proxy for total linear polarization. (Ru) is 
only used to help correct for errors in aircraft geometry. 

Our radiative transfer model uses Mie theory to compute 
the polarized phase functions (single scattering) for vari- 
ous types of aerosols or cloud droplets, and Rayleigh scat- 
tering for gases (for a historical review of Mie theory, see 
Mishchenko and Travis, 2008). This means we require a 
measure of cloud droplet size distribution, aerosol size distri- 
bution (which is typically bimodal) and the complex aerosol 


refractive index ((m — (m) + S5(m))), and that we assume 

the aerosols scatter like spheres. We express size distribution 
as an effective radius, (r e ), and variance, (v e ), where the for- 
mer has units of microns and the latter is unitless (Hansen 
and Travis, 1974). Integrals of the size distribution are per- 
formed within a 0 to 20 pm range. The size distribution of 
aerosols is typically bimodal, so we use the subscripts (/) 
and (c) to refer to fine and coarse mode aerosols. The for- 
mer is usually the product of a chemical reaction, such as 
sulfate photo-oxidation or the production of biomass burning 
aerosols, and the latter is due to physical processes such as 
the creation of sea salt or dust aerosols. Cloud size parame- 
ters are denoted with the subscript (cl). We use the lognor- 
mal distribution for aerosols and the gamma distribution for 
cloud droplets (equations 2.60 and 2.56 in Hansen and Travis 
(1974)). The complex refractive index, (m), is set indepen- 
dently for each mode. In this work we assume the real part 
of (m) is spectrally independent, while the imaginary com- 
ponent of the fine mode is specified with two parameters (see 
Sect. 2.4). 

Once the single scattering properties for an aerosol or 
cloud droplet have been determined, multiple scattering is 
computed for the full atmosphere to create the polarized re- 
flectance at the observation altitude. We use the Doubling 
and Adding technique (Hansen and Travis, 1974; De Haan 
et al., 1987), which assumes plane parallel layers of clouds 
and aerosols. The altitude of these layers are determined dur- 
ing the J-31 atmospheric profile as described in Sect. 2.2 
and they are held constant in the optimization. The op- 
tical depth, (r), of each layer is also required. As noted 
above, initial aerosol optical depth values are derived from 
the AATS-14. Cloud optical depth is set arbitrarily large 
((r c i(0.555 pm) = 10.0)) to saturate the signal, since polarized 
reflectance is only sensitive to the first few units of cloud op- 
tical depth. Because of this, the surface is obscured and we 
assume its reflectance is negligible. 

A complete list of the model values is presented in Table 1. 
Italicized parameters are allowed to vary during optimization 
and are thus retrieved, while other parameters are fixed and 
considered part of the radiative transfer model. 

2.1 AAC scene 

On 13 March 2006, the J-31 aircraft, containing the RSP and 
a variety of other instruments, departed from the Veracruz 
airport. The airplane flew northwest along the Gulf of Mex- 
ico coast and above a marine stratocumulus cloud. The cloud 
was in the process of dissipating, and it had a larger geo- 
graphical extent earlier in the day (Bergstrom et al., 2010). 
Following the cloud overflight, the J-31 descended to the sur- 
face in a recently cloud free region. Figure 2 is a map of 
this scene, the descending spiral, and the surrounding area. 
We used observations that were made in that spiral profile 
to constrain cloud and aerosol layer heights in our scene, on 
the assumption that the atmospheric vertical profile had not 
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Table 1. Retrieved aerosol parameters and the initial values used in 
optimization. Parameters listed in italics are optimized during the 
retrieval, while all others are kept constant. * : Refractive Index val- 
ues have no spectral dependence in the initial values, although the 
spectral dependence of the fine mode aerosol imaginary component 
is allowed to vary as described in Sect. 2.4. 


Parameter 

Symbol 

Initial value 

Fine mode aerosols 



Optical depth, 0.555 pm 

T f (0.555) 

0.12 

Real refractive index* 

1 R(m f) 

1.47 

Imaginary ref. index, 0.555pm* 

3(m f) 

0.014 

Imaginary refractive index, p* 


0.007 

Imaginary refractive index, q* 

%(mf, q ) 

0.007 

Effective radius, (pm) 

r e,f 

0.20 

Effective variance 

v e,f 

0.20 

Coarse mode aerosols 

Optical depth, 0.555 pm 

r c (0.555) 

0.04 

Real refractive index (*) 

rnm c ) 

1.47 

Imaginary refractive index (*) 

k c 

0.014 

Effective radius (pm) 

r e.c 

7.67 

Effective variance 

Ve.c 

0.49 

Cloud droplets 

Optical depth, 0.55 pm 

t c1 (0.555) 

10.0 

Effective radius (pm) 

r e,cl 

6.25 

Effective variance 

v e,c\ 

0.075 


changed dramatically in the 125 km distance between the ob- 
servation location and descending spiral. Figure 3 is the tem- 
perature profile and vertical change in aerosol optical depth 
from this spiral. The temperature profile was used to deter- 
mine the cloud top height (480 m) for our radiative transfer 
model, as shown in blue in Fig. 3. Change in aerosol optical 
depth (equivalent to the aerosol extinction coefficient), from 
the AATS-14 instrument, was used to constrain the AAC 
height (550 m to 750 m), shown in grey. Table 2 lists the 
geometry and other specifics associated with this scene and 
the downward spiral. A full description of this and other J- 
31 flights during MILAGRO can be found in Molina et al. 
( 2010 ). 

2.2 Aerosol sources 

The AAC we encountered in our scene were a complicated 
mixture with a variety of sources. Wind conditions on the 
13th of March were light and variable, preceded by south- 
westerly winds associated with a trough over the western 
portion of the USA. This portion of the MILAGRO field 
campaign was characterized by low humidity and few cir- 
rus clouds, and there were probable interactions between 
MCMA emissions and biomass burning smoke (Fast et al., 
2007). Aerosols from a variety of sources were found dur- 
ing MILAGRO, including sulfates of both industrial and vol- 


Table 2. Geometry and other parameters associated with our AAC 
scene, and the subsequent downward spiral. Aerosol optical depths 
are from the AATS instrument 


Parameter 

Value 

AAC scene 

Starting Time 

16:05 UTC 

J-31 altitude 

5380 m 

Latitude 

20.15° N 

Longitude 

96.68° W 

Solar Zenith Angle 

44° 

Relative Azimuth Angle 

20° 

Number of averaged scans 

22 

Downward Spiral 

Starting Time 

16:22 UTC 

J-31 altitude 

5360 to 51 m 

Latitude 

21.06° N 

Longitude 

96.95° W 

Cloud top 

480 m 

Aerosol layer bottom 

550 m 

Aerosol layer top 

750 m 

(r0.519 pm) above 55 0m 

0.165 

(r0.604pm) above 550 m 

0.130 

(t 2.139 pm) above 550m 

0.040 


canic origin (de Foy et al., 2009; Grutter et al., 2008), and 
organic carbon. Secondary organic carbon was more domi- 
nant than primary organic carbon, and was the result of both 
urban and (generally anthropogenic) biomass burning activi- 
ties (Doran et al., 2007; Crounse et al., 2009; Yu et al., 2009; 
Paredes-Miranda et al., 2009). The MCMA is in a high valley 
(2200 m above sea level), so it is feasible for anthropogenic 
aerosols to be lofted above low level marine stratocumulus 
clouds on the coast. 

This leads us to believe that the AAC aerosols in our scene 
are probably anthropogenic and that the fine size mode is 
the dominant contributor to scattering and absorption. To 
further investigate the aerosol sources, we used the FLEX- 
PART Lagrangian particle dispersion model (Stohl et al., 
2005) to determine the sensitivity of our scene to emission 
in other regions. Figure 4 is the footprint emission sensi- 
tivity obtained from FLEXPART driven by European Cen- 
tre for Medium Range Weather Forecasts (ECMWF) mete- 
orological data. Values of high emission sensitivity can be 
found both over the MCMA and biomass burning areas in 
the vicinity of Veracruz. However, results are somewhat dif- 
ferent when the meteorology is driven by an alternate model, 
the Global Forecast System (GFS) from the National Center 
for Environmental Prediction (NCEP). FLEXPART driven 
by GFS shows an increased influence of Veracruz area smoke 
at the expense of MCMA pollution. This disparity is prob- 
ably due to the weak and variable winds at that time, and 
indicates that the aerosols have had the opportunity to mix. 
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MODIS Terra image from 17:20 UTC 


Gulf of Mexico 


RSP scene 

16:05 UTC; Altitude: 5380m 
Relative Azimuth angle: 20° 
Solar Zenith Angle: 44° 
20.15° N. 96.68° W 


AERONET sites 


Descending spiral 


J-31 Aircraft flight track 
(South to North) 


Mexico City 


MODIS fire locations 
(previous eight days) 


Fig. 2. The spatial context of our AAC scene is presented in this figure. The blue circle indicates the location of RSP observations, above a 
marine stratocumulus cloud on the Gulf of Mexico coast. A portion of the J-31 flight track is shown in yellow. The J-31 performed a spiral to 
the surface about 125knr northwest of the scene, and data collected during this descent provided information about cloud and aerosol vertical 
distribution. Aerosol sources include urban/industrial emissions in the Mexico City Metropolitan Area (MCMA) basin, a high valley to the 
west, and numerous (mostly agricultural) fires indicated by the red fire icons. Fire locations were identified by the MODIS active fire product 
and represent fires within the previous eight days. The MODIS Terra instrument captured the underlying image about an hour and fifteen 
minutes after our scene was observed. 


accumulate, and age, with the exact distribution of contribut- 
ing sources being somewhat more uncertain than for other 
cases. 

2.3 Determination of initial cloud size distribution 


Prior to performing our retrieval, we made an initial estimate 
of the cloud size distribution using a Look Up Table (LUT) of 
cloud droplet single scattering properties. Goloub, in his pa- 
per describing cloud property retrieval from the Polarization 
and Directionality of Earth Reflectances (POLDER) instru- 
ment (Goloub et al., 2000), found that polarization of cloud 
reflectances is less sensitive to multiple scattering than total 
reflectance. As cloud optical thickness becomes greater than 
2, ( Rq ) in the single scattering approximation becomes 




1 

4(cos0^ +cos0„) 


Q{^> r e,cl> fe,cl) 


( 2 ) 


where ( 9 S ) is the solar zenith angle, ( 6 V ) is the view zenith 
angle and ( Q ) is the Stokes vector component of the single 
scattering phase function. (£) is the scattering angle, defined 
with respect to the incident illumination vector, and (r e . c |. 
u e , c i) are the cloud effective radius and variance. 


At large scattering angles, where (Rq) is dominated by 
the cloud signal (see Fig. 1), we can identify the cloud size 
distribution independently from cloud optical thickness or 
multiple scattering effects. (Rq) is sensitive to cloud effec- 
tive radius at scattering angles larger than 130°, and effective 
variance at angles larger than 145° (Waquet et al., 2009b). 

Our initial estimate of cloud size distribution is determined 
by matching (Rq) to a LUT of ( g(£, r c c i, fe,cl)) for a vari- 
ety of size distributions. We used scattering angles between 
135° and 162° (which was the largest scattering angle avail- 
able in our scene). The LUT contained effective radii be- 
tween 5 and 20 pm at 0.25 intervals, and effective variances 
between 0.01 and 0.25. Aerosol and molecular effects above 
the cloud are expressed in a manner that does not change 
rapidly with scattering angle (unlike the sharp cloud bow fea- 
tures), so we used a Fourier domain high pass filter (Gonza- 
lez and Woods, 1992) to remove this low frequency signal 
prior to matching (Rq) to the LUT. To account for the possi- 
bility that absorbing AAC dampen the cloud signal, we allow 
both (r e , c i) and (u e , c i) to vary during optimization. 
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Temperature (C) 

20 21 22 23 24 25 



Aerosol Extinction at 519nm [1/km] 

Fig. 3. The aerosol optical depth and temperature during the down- 
ward spiral northwest of our scene are presented in this figure. Opti- 
cal depth observations were made by the AATS-14 instrument, and 
the derivative with respect to height (equivalent to the extinction 
coefficient) of observations at 0.519 pm is presented in black. Two 
layers of aerosols are present. The lowest, between about 150 and 
350 m, was most likely not observable in our scene. The upper layer, 
between 550 and 750 m, was probably above the cloud and therefore 
the same AAC that we observed in our scene. The atmospheric tem- 
perature profile is shown in green. At 48 m, there is a sharp change 
in the temperature profile. This probably marks the upper bound- 
ary of the marine stratocumulus cloud that existed at this location 
several hours before. Based on this information, we constructed our 
radiative transfer model so that the cloud top is at 480 m, with an 
aerosol layer above between 550 and 750 m. The initial aerosol op- 
tical thickness value of the upper layer was estimated to be 0.16 at 
0.555 pm from AATS-14 data. 

2.4 Optimization 

Aerosol and cloud optical properties are retrieved by com- 
paring a radiative transfer simulation of our scene to RSP 
observations, and tuning the model parameters until the best 
match is found. We use the Levenberg-Marquardt optimiza- 
tion technique as described in More (1977) and Markwardt 
(2009), and implemented for the RSP in Knobelspiesse et al. 
(2011). This means the retrieved parameters are not real 
numbers but regions of state space that are associated with 
the scene considering observational uncertainty. Waquet 
et al. (2009a) also used a similar approach with RSP data 
for aerosol retrievals over land. 


Footprint emission sensitivity in the nested domain for 2006/03/13 

Lower release height: 500m; Upper release height: 750m 
Meteorological data are from ECMWF 


Fig. 4. Emission sensitivity footprint for the 500-750 m altitude 
range at our AAC scene obtained from FLEXPART driven by me- 
teorological data from ECMWF. 

We used the software developed in Knobelspiesse et al. 
(2011) for this study with slight modifications. First, we 
increased the numerical accuracy of the radiative transfer 
model. While computationally more expensive, this is re- 
quired to properly simulate cloud droplets, which have much 
stronger forward scattering and produce more dramatic fea- 
tures than aerosols (such as cloud bows). We also modified 
the software so that the aerosol optical depth of a layer is 
a directly retrieved parameter. Previously, the total quantity 
of aerosols in a layer was determined by the particle number 
concentration, and this was the parameter that was tuned dur- 
ing optimization. Validation, however, is usually performed 
by comparisons of optical depth, since this is a common 
observation made by instruments such as sun photometers. 
While optical depth can of course be computed using the 
number concentration (and aerosol extinction cross section, 
which depends on size and refractive index), we switched 
to directly optimizing this parameter so that error computa- 
tion is simpler. We also noticed that the optimization itself 
appears to require fewer iterations and is better behaved, al- 
though this observation is limited to this AAC scene. 

We also modified the software in Knobelspiesse et al. 
(2011) to allow the imaginary component of the refractive 
index to have a spectral dependence. Rather than using the 
single value for each size mode, we now use a pair of param- 
eters, (%(m p )) and (S(m 9 )): 

»(m(A.)) = »(OT p ) + S(m 9 )(k + 0.445)“ 10 (3) 

where (X) is the wavelength in microns. This spectral de- 
pendence function is intended to allow for the increase in 



0.000 0.001 0.002 0.005 0.010 0.020 0.039 0.078 0.156 0.313 0.625 1.250 

ns / kg 

Maximum value: 8.22 ns/kg 
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absorption at shorter wavelengths due to “brown carbon” 
(BrC, absorbing organic carbon aerosols). This was observed 
by the SSFR instrument on the same flight as our AAC scene 
(Bergstrom et al., 2010), which is our primary motivation for 
this change. 

The choice of this form of specular dependence is unfortu- 
nately somewhat arbitrary, as direct observations of the imag- 
inary refractive index of aerosols in their entirety are difficult 
and quite limited. Evidence of absorption in observations of 
the aerosol optical depth, on the other hand, are much more 
common. The absorption optical depth, (r a ), which is the op- 
tical depth due to absorption (but not scattering) is generally 
characterized with the equation 


Imaginary refractive index spectral dependence 

RSP fit resuits p.O 0000 q: 0.0031 

Rowers et al. (2010) Low OC p:0.0068 q: 0.0010 
Flowers et al. 2010) Medium OC p:0.01 1 1 q: 0.0012 
Flowers et al. (2010) High OC p:0.0141 q: 0.0029 

Sun et al. (2007) OC1 model p:0.0007 q: 0.0045 



0.6 0.7 0.8 

Wavelength, microns 


t«W = r«,»(W“ (4) 

where (X 0 ) is a reference wavelength and ( r a o ) is the ab- 
sorption aerosol optical depth at that wavelength, (a) is the 
absorption Angstrom Exponent, which should be equal to 1 .0 
for very small particles with spectrally invariant refractive in- 
dices (vandeHulst, 1981 and Bond, 2001). Larger Angstrom 
Exponent values have been found in some biomass burning 
(smoke) aerosol plumes, such as in Lewis et al. (2008), indi- 
cating an increase in the imaginary refractive index (and thus 
absorption) at shorter wavelengths. Authors such as Hoffer 
et al. (2006) and Dinar et al. (2008) have found imaginary 
refractive indices that increase with decreasing wavelength 
for Humic -Like Substances (HULIS), which are likely com- 
ponents of BrC, but direct observations of entire aerosol par- 
ticles are limited. Others (Chen and Bond, 2010, Kirchstet- 
ter et al., 2004 and references therein) have found increas- 
ing aerosol absorption coefficients with shorter wavelength. 
While this is evidence of an increase in shortwave imaginary 
refractive index, the absorption coefficient is also dependent 
on particle size and real refractive index. 

For small particles, the imaginary refractive index 
is proportional to wavelength and absorption coefficient 
((2(m(k)) ocka)) (Sun et al., 2007), which we use to con- 
nect our model for imaginary refractive index in equation 
3 to observations. The (3(w p )) parameter indicates the 
quantity of spectrally independent absorption, while (3 (m q )) 
expresses the quantity of spectrally dependent absorption. 
A large negative exponent has been selected so that varia- 
tions in and (S(m 9 )) have very different impacts, 

which allows their simultaneous retrieval during optimiza- 
tion, since they are more orthogonal in retrieval space. The 
exponent value of (—10) was chosen from the upper limit 
to the absorption A ngstrom Exponent found by Chen and 
Bond (2010). Spectral dependence is also shifted so that 
(k(0.555) = £5(m p ) + 2s(m ? )) and converges to 

(2(m p )) to as wavelength increases. 

Figure 5 presents the imaginary refractive index that was 
retrieved in this work (in black, see Sect. 4), along with the 
imaginary refractive index models representing other scenes 
from Flowers et al. (2010) and Sun et al. (2007) (presented 


Fig. 5. Imaginary refractive index spectral dependence models from 
Flowers et al. (2010) (red, green and blue) and Sun et al. (2007) 
(magenta), compared to RSP retrieval results (black, see Sect. 4). 
Diamonds represent the published models, while dashed lines are 
the fit from equation 3 to that model. 

as diamonds). Equation (3) was fit to these models, and is 
plotted as dashed lines in Fig. 5. Our imaginary refractive in- 
dex parameterization is clearly capable of representing sev- 
eral types of refractive index models, and is therefore suitable 
for use during optimization. 

Our optimization methods provide the retrieval error for 
each model parameter (see Sect. 3). Since it is not directly 
parameterized, the uncertainty in ( k ) must be derived from 
the uncertainties in (%(m p )) and (£y(m 9 )). Assuming that 
correlation between (%(jn p )) and is minimal, (o>) 

for those parameters 

o* 2 W = ol p + <*l q (k + 0.445)- 20 . (5) 

3 Simulated aerosol above cloud 

To help us design an appropriate optimization strategy, we 
first assessed the information content available in an AAC 
scene. We did this by simulating a climatologically defined 
aerosol at a variety of optical depths suspended above a ma- 
rine stratocumulus cloud. Other than optical depth, aerosol 
optical properties were chosen to match the “Mexico City Ur- 
ban/Industrial” class described in Dubovik et al. (2002) and 
used as the initial value during optimization. 

For each simulation, we numerically estimated the Jaco- 
bian matrix, (J), which expresses the sensitivity of the for- 
ward model to change in each parameter. 


where the artial derivative of forward model, (F), for the sim- 
ulated set of parameters, (x), is computed for each observa- 
tion ((/)) and each parameter ((/)). The partial derivative 
was estimated numerically by perturbing the ( j ,h ) element of 
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(x) (which we denote as (x'j) and recalculating the forward 
model, (F). 


JiM)' 


Fi&)-Fi(x) 


x j ~ X J 


(7) 


While the Jacobian is a useful metric to identify the rel- 
ative importance of different parameters in measurement 
space, we are interested in parameter sensitivity. This is ex- 
pressed with the retrieval error covariance matrix. 


C, = (J r Cr 1 J)" 1 


( 8 ) 


This requires the measurement error covariance matrix, 
(C7-), which is specific to the RSP and accounts for measure- 
ment uncertainty due to errors in calibration, and observation 
geometry. It is computed as in Eq. (3) of Knobelspiesse et al. 
(2011). The square root of the diagonal elements of (C v ) are 
the standard deviations of the errors for each parameter in 
(x), provided that (C t) is accurate and the forward model is 
linear over the perturbation range used to numerically calcu- 
late (J) (Hasekamp and Landgraf, 2007). 

We simulated a scene very similar to the initial values we 
selected for optimization, with the exception of optical depth. 
Aerosol and cloud optical properties that were used for the 
simulation are listed in Table 1, while the cloud and aerosol 
vertical distribution is described in Table 2. Unlike our RSP 
scene, however, we selected a view zenith angle of (45°) and 
relative azimuth angle of (45°). This geometry was chosen to 
generalize the results for the APS, which will typically make 
observations farther from the solar principal plane (where 
(ip = 0°)) than the RSP did in this scene. Observations in the 
solar principal plane contain more information than those at 
other relative azimuth angles because a larger range of scat- 
tering angles are measured. For this reason, we would expect 
lower errors for the RSP viewing geometry, although errors 
from brief simulations we performed at observational relative 
azimuth angles were nearly identical. 

We performed two different sets of simulations. The first 
set simulated an optimization that used (Rq) in seven RSP 
bands (centered at 0.410, 0.470, 0.555, 0.670, 0.865, 1.590 
and 2.250 pm) at half degree intervals in the view zenith an- 
gle range of (±60°) from nadir. At this geometry, the scatter- 
ing angle range is from (85°) to (150°) and includes both side 
scattering angles (where aerosol properties dominate the sig- 
nal) and the cloud bow (where cloud properties dominate the 
signal). The second simulation also used (Rq) at the wave- 
lengths described above, but view angles restricted to (20°) 
to (60°) in the forward direction. The scattering angles that 
correspond with this are between (85°) and (120°), excluding 
the cloud bow. This is similar to the scattering angle range 
used to retrieve AAC properties in Waquet et al. (2009b), 
who was limited in his ability to determine cloud proper- 
ties (due to the coarser angular resolution of POLDER) and 
performed an optimization using only observations where 
aerosol properties dominate. 


Simulated retrieval errors are shown in Fig. 6. Aerosol 
remote sensing accuracy requirements for climate models, 
as assessed by Mishchenko et al. (2004), are also presented. 
These requirements were not determined for AAC, but for re- 
trievals of aerosols over land (without the presence of clouds) 
or clouds without aerosols above. We include them in Fig. 6 
to compare the potential of AAC retrievals with the stan- 
dard retrievals by the APS instrument for which these ac- 
curacy requirements were determined. For the full angu- 
lar range simulation, these accuracy requirements are indeed 
met for aerosol optical depth, fine mode aerosol size, and 
cloud droplet size. Real refractive index and single scat- 
tering albedo requirements are not met, although errors de- 
crease and nearly approach the requirements as aerosol op- 
tical depth increases. Errors for the simulation with fewer 
viewing angles are clearly higher. Aerosol optical depth ac- 
curacy requirements are no longer met, and errors for fine 
mode effective variance are too high at very low optical 
depths. Simulations also show sensitivity to cloud droplet 
size and effective variance within accuracy requirements for 
aerosol optical depths less than one. This is surprising 
considering that this simulation did not use observations of 
the cloud bow, where cloud properties should be dominant. 
These results, however, do not indicate the degree of corre- 
lation between different parameters in the retrieval. For the 
simulation with a total optical depth of (r (0.555) = 0. 12), the 
cloud effective radius was strongly correlated with the fine 
mode optical depth and fine mode effective variance (see ap- 
pendix A). Cloud effective variance was strongly correlated 
with the optical depth of both fine and coarse aerosols. These 
strong correlations do not exist for the equivalent retrieval 
that used the full angular range of observations. 

These simulations suggest several strategies for success- 
ful optimization. The most obvious is that optimizations 
should utilize observations both of the cloud bow and side 
scattering angles. Cloud size parameters should be retrieved 
during optimization, since their assessed errors are low and 
those parameters are not significantly correlated with oth- 
ers (with one exception, see Appendix A). Furthermore, fine 
and coarse mode aerosol optical depth and fine mode size 
parameters can all be retrieved accurately. The real refrac- 
tive index and absorption, however, are not retrieved accu- 
rately. Furthermore, those parameters have a high degree 
of correlation between themselves and with other parame- 
ters (see Appendix A). An ideal situation would be to have 
an accurate external measure of the aerosol complex refrac- 
tive index that can be specified during optimization. Obvi- 
ously, such measurements typically do not exist, so we are 
forced to consider retrieving some of the complex refractive 
index parameters, even if their retrieval error is larger than 
requirements for typical aerosol values. We are most inter- 
ested in AAC absorption, since absorption is an important 
factor in the overall radiative forcing of AAC. The imagi- 
nary component of refractive index, however, is highly cor- 
related with the optical depth of coarse mode aerosols. For 
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Fine Mode Optical Depth, 555nm 


Fine Mode real refractive index 


Fine Mode effective radius 






Coarse Mode Optical Depth, 555nm 





AOT(555nm) 


Fine Single Scattering Albedo, 532nm 


o.i 1.0 

AOT(555nm) 



Dashed: MC, side scattering only 
Dotted: error thresholds 


Fig. 6. Simulated retrieval error for “Mexico City” type aerosols (see Dubovik et al., 2002) suspended above a marine stratocumulus cloud. 
The aerosols, which were placed in a uniform layer between 550 and 750 m above sea level, were simulated at a variety of optical depths. 
The x-axis of each plot is the fine mode aerosol optical depth at 0.555 pm. Coarse mode optical depth was maintained at 20 % of the fine 
mode value for each simulation. The marine stratocumulus cloud had an effective radius of (r e c ] = 6.25 pm) and effective variance of 
(n e c ] = 0.075), an optical depth of 10, and was uniformly distributed (with respect to pressure) from the ground to 480 m. Solid lines are 
the simulated errors for an observation including all view zenith angles within (60°) of nadir, which is similar to our scene and includes the 
cloud bow (and corresponds to a scattering angle range of (85°) to (150°)). Dashed lines are the simulated errors for an optimization that 
only uses observations at angles between (20°) and (60°) in the forward direction, which excludes the cloud bow and is similar to the method 
of Waquet et al. (2009b), who had greater uncertainty when determining cloud parameters. These angles correspond to a scattering angle 
range of (85°) to (120°). Dotted lines are the accuracy thresholds described in Mishchenko et al. (2004) required to improve aerosol climate 
models. It is important to note that these values are intended for retrievals of aerosols in the absence of clouds or clouds in the absence of 
aerosols - not AAC. Also, the threshold for effective radius (top right) is (0.1 pm), greater than all the simulated errors in that plot. The single 
scattering albedo (bottom, center) accuracy threshold is 0.03, less than all simulated error values. 


our optimization, we decided to fix the real refractive index 
at the climatological value for Mexico City aerosols as con- 
tained in Dubovik et al. (2002), while allowing both parame- 
ters describing the imaginary refractive index (see Sect. 2.4) 
to vary. Since the latter are correlated with coarse mode 
optical depth, we fixed that parameter to the optical depth 
observed by the AATS sun photometer at long wavelengths 
(where the coarse mode should contribute the most to overall 
extinction). Free parameters in our retrieval are therefore the 
aerosol fine mode optical depth, both parameters describing 


the imaginary component of fine mode refractive index, fine 
mode effective radius and variance, and cloud droplet effec- 
tive radius and variance. 

4 Results 

4.1 Cloud properties 

As described in Sect. 2.3, initial values for cloud droplet size 
were determined by matching a look up table of cloud droplet 
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Cloud fit to observation 



Fig. 7. RSP observations (solid lines) were matched to a look up 
table of cloud optical properties (dashed lines) to determine ini- 
tial cloud droplet size properties for the optimization. The cloud 
that matched best had droplets with an effective radius of (r e c [ = 
6.25 pm) and effective variance of (u e c i = 0.075). Low frequency 
differences (bias) between cloud and RSP ( Rq ) are due to the ef- 
fects of aerosols and Rayleigh scattering on the observation. 


single scattering properties to observations. Because RSP 
observations also contain the effects of Rayleigh scattering 
and aerosols above the cloud, the cloud parameters deter- 
mined with this method were allowed to vary in the next 
step, optimization. Here, aerosol and Rayleigh scattering ef- 
fects were minimized by removing the low angular frequency 
trend before comparing the observation to the look up table. 
The best match was found for a cloud whose droplets have 
an effective radius of (r e , c i = 6.25 pm) and effective variance 
of (u e ,ci = 0.075). Fig. 7 shows this match (before low fre- 
quency trends are removed) for three wavelengths. While 
there are significant offsets between the RSP observation and 
cloud model, the angular location of the cloud bow in the data 
are well represented with the model. 

4.2 Aerosol and cloud simultaneous retrieval 

Optimization was performed using the data, methodology 
and initial values described above to retrieve seven aerosol 
and cloud parameters. Eleven iterations of the Levenberg- 
Marquardt method were required to converge to a solution. 
Fig. 8 shows the radiative transfer model match to the ob- 
servations, while table 3 is the aerosol and cloud parameters 
that were retrieved from the model. Comparisons between 
these results and observations made by other instruments are 
presented in Fig. 9. 


AAC optimization results 



Fig. 8. RSP observations (solid lines), and the model match (dashed 
color/black lines) to those observations. 


The optimized radiative transfer model appears to match 
best at forward scattering angles, where aerosol properties 
dominate. This is somewhat understandable in that the mea- 
surement error (as expressed in the error covariance matrix, 
(C 7 -)) is smallest for these viewing angles. Measurement 
error is largest for cloud bow view zenith angles (between 
(—20°) and (0°)) because of uncertainties in aircraft geom- 
etry that have the largest impact on features that change 
rapidly with angle. Perhaps as a result, model results are 
worst here, and underestimate the polarization in the cloud 
bow in all channels. However, we also found that few com- 
binations of aerosol or cloud parameters could create such 
strong polarization in the cloud bow, indicating that our 
cloud model may be too simplistic. Regardless, aerosol and 
cloud parameter values have errors that are similar to those 
of the simulation in Sect. 3 for comparable optical depths. 
Aerosol optical depth and size parameters have low errors 
that meet the accuracy requirements of Mishchenko et al. 
(2004). There are no accuracy requirements for the aerosol 
imaginary refractive index, but the related single scattering 
albedo has an uncertainty so large as to render that parame- 
ter nearly useless. This is not much of a surprise, however, 
since the predicted error in single scattering albedo in the 
sensitivity test is comparable. The sensitivity test also shows 
that this error will decrease with increasing aerosol optical 
thickness, so that optical thicknesses above (r (0.555) =0.8) 
have single scattering albedo errors nearly within accuracy 
requirements. The global radiative impact of these errors will 
therefore depend on the AAC optical thickness distribution. 

The retrieved AAC properties show a somewhat thin 
(about half the climatological median optical thickness) layer 
of aerosols that have a fine mode effective radius very simi- 
lar to the climatology of Dubovik et al. (2002). This is rea- 
sonable, however, because the observed altitudes represent 
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Table 3. Retrieved aerosol and cloud parameters. *: These pa- 
rameters are derived from other parameters, not directly retrieved. 
Errors in bold are larger than accuracy guidelines from Mishchenko 
et al. (2004) (note that guidelines were not provided in this publica- 
tion for imaginary refractive index, and were not intended for AAC 
retrievals). 


Symbol 

Unit 

Retrieval 

Error 

Simulation Error 

ry (0.555) 

n/a 

0.104 

0.016 

0.009 

r(0.555)* 

n/a 

0.144 

0.016 

0.009 

S( mf' P ) 

n/a 

0.000 

0.062 

0.064 


n/a 

0.003 

0.036 

0.015 

S(m /X0.532)* 

n/a 

0.004 

0.077 

0.063 

r e,f 

(im 

0.143 

0.007 

0.029 

v e,f 

n/a 

0.057 

0.036 

0.106 

r e,c\ 

(jm 

6.82 

0.187 

0.061 

v e,c\ 

n/a 

0.028 

0.009 

0.012 

m-(0.532)* 

n/a 

0.868 

0.450 

0.306 


only a portion of the atmospheric column and do not in- 
clude, for example, boundary layer aerosols. The width of 
the fine mode size distribution is quite small, which may be 
explained by the hypothesis in Bergstrom et al. (2010) that 
these aerosols have been cloud processed (although it should 
be noted that aerosol size distribution was not directly ob- 
served in that paper). The cloud droplet effective radius is 
close to, but slightly smaller than, the climatology for ma- 
rine clouds in Miles et al. (2000). Compared to the initial 
values, the cloud effective radius increased slightly, while the 
effective variance decreased. This is somewhat understand- 
able in light of the method used to determine the cloud initial 
size parameters. Aerosols and Rayleigh scattering above a 
cloud would generally be expected to attenuate the strong 
polarized signal emanating from the cloud bow. This would 
act to decrease the polarization strength of the cloud bow, 
and dampen the secondary oscillations associated with nar- 
row size distributions. Once incorporated into a full multiple 
scattering radiative transfer model, the effects of droplet size 
on the observations are more accurately calculated, provid- 
ing a better estimate of the droplet size distribution. 

One interesting result from this optimization is that it con- 
verged to a solution where the (%(m f q )) parameter is dom- 
inant, while Q(m f p )) became nearly zero. Initial values 
for both these parameters were an equal split of (£s(m/)) 
from the Dubovik et al. (2002) Mexico City climatology. 
Errors for both (3(w f p )) and (3(w f q )) are large, but this 
indicates a likely spectral dependence in the imaginary re- 
fractive index. This was also found by Bergstrom et al. 
(2010) in a nearby scene. We also tested optimizations that 
had spectrally flat initial conditions ((£s(m f, p ) — 0.014) and 
(3(m f q ) = 0.0)), but they converged to what we found here, 
only with many more Levenberg-Marquardt iterations. 


Spectral Aerosol Optical Thickness 




Fig. 9. The topmost panel in this figure is the spectral dependence of 
retrieved total aerosol optical thickness (black line) and the AATS 
observation (green line) at an altitude of 480 m during the down- 
ward spiral indicated in Fig. 2. The bottom panel is the spectral 
dependence of the total single scattering albedo. The black line is 
the RSP retrieved single scattering albedo, where the large error 
bars are omitted for clarity (they would fill the entire range of val- 
ues in the plot). Red indicates SSFR observations as described in 
Bergstrom et al. (2010). 

4.3 Comparison with other observations 

Observations of the aerosol plume were also made by other 
instruments on the nl-31 aircraft. The Ames Aerosol Track- 
ing Sunphotometer (AATS) made observations of the spec- 
tral aerosol optical depth in a cloud free region about 125km 
northwest of our scene. Aerosol absorption was also ob- 
served in this area by the Solar Spectral Flux Radiometer 
(SSFR) (Pilewskie et al., 2003; Schmidt et al., 2010), and 
published in Bergstrom et al. (2010) (see Fig. 2). As de- 
scribed in Sect. 2.2, atmospheric transport was weak prior to 
our observation, so we assume the aerosols at this location 
are similar to those in our AAC scene. 

The aerosol optical thickness we retrieved agrees very well 
with AATS observations from an altitude of 480 m (which 
was the cloud top height in our AAC scene). To some ex- 
tent, this is by design, since the longest wavelength AATS 
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optical thickness values were used to constrain the AAC 
coarse mode aerosol optical thickness during optimization. 
Even so, at wavelengths less than (0.8 |jm), the majority of 
total aerosol optical thickness is determined by the fine mode, 
and AATS measurements at these wavelengths agree well 
within retrieval uncertainties. 

The single scattering albedo that we retrieved has a very 
large uncertainty. However, comparisons with SSFR ob- 
servations show some degree of similarity, at least for the 
shorter wavelengths of the spectrum. It is encouraging that 
agreement is better at these wavelengths since this is where 
the fine aerosol mode dominates extinction, and the imag- 
inary refractive index of the fine mode was allowed to vary 
during optimization. Differences are greatest for the (0.7 pm) 
to (1.1 pm) wavelength range. The coarse aerosol mode con- 
tributes more to the total extinction at these wavelengths, 
but the imaginary refractive index for the coarse mode was 
held fixed during optimization. This may indicate that the 
coarse mode imaginary refractive index was not properly se- 
lected. Given the overall imaginary refractive index uncer- 
tainty, however, no firm conclusions can be made. 

4.4 Model uncertainty 

Forward model suitability is a major issue for optimizations 
such as this. Because of the finite information content avail- 
able in our observations, we were required to make several 
assumptions about the nature of our scene. Many of these as- 
sumptions were based on external observations (such as the 
aerosol vertical distribution from AATS) or climatologies of 
aerosol properties (such as the coarse mode aerosol optical 
properties from Dubovik et al., 2002). While it is difficult 
to quantitatively assess these assumptions without further in- 
formation, we can determine if the uncertainty in the data 
that went into making them has an impact on the retrieved 
parameter values. To do so, we start by modifying one of 
the assumptions in the forward model by its uncertainty, and 
recomputing with the retrieved parameters. The difference 
between this (which we will denote (F(x))) and the forward 
model results is the model error in measurement space. We 
can then use the Jacobian to project the error ((e)) into state 
space. 

e=(JJ r r 1 J[F(x)-F(x)] (9) 

Table 4 contains values of (e) as a percentage of the re- 
trieved parameter value for various model assumptions. 

Model uncertainties come from a variety of sources. One 
class of model assumptions have to do with the physical de- 
scription of the aerosol and cloud scene as plane parallel lay- 
ers in a radiative transfer model. In this paper we test the 
validity of the vertical layering heights. Validity of the plane 
parallel model approach requires a comparison with a three 
dimensional radiative transfer model, which is beyond the 
scope of the current work and will be a topic of future re- 
search. Another class of assumptions has to do with the ge- 


ometry of the scene. While presumably this will not be an 
issue for orbital instruments, the aircraft carrying the RSP 
is subject to rapid attitude changes, so the measurement ge- 
ometry does have uncertainty. Finally, we also test the as- 
sumptions we make about the aerosol optical properties we 
were not able to retrieve, such as the fine mode real refractive 
index and the coarse mode size and refractive index. 

Aerosol and cloud layer tops were perturbed by the ver- 
tical resolution of the CAFIPSO lidar (Winker et al., 2006 
and Winker et al., 2003). CAFIPSO is in the “A-Train” or- 
bit, which was also the destination of the APS instrument. 
Operational retrievals of AAC from instruments similar to 
APS would likely use the cloud top height and aerosol ver- 
tical distribution from CAFIPSO. Therefore, we used 30 m 
perturbations, which are equivalent to the vertical resolu- 
tion of CAFIPSO for the lower troposphere. It is clear that 
(e) for most retrieved parameters is not strongly affected by 
this level of uncertainty in the vertical distribution, although 
cloud layer top is more important than the aerosol vertical 
distribution. The imaginary refractive index for the fine mode 
does have a somewhat large degree of error. However, this er- 
ror is much smaller than the retrieval error itself, which was 
about twenty times larger than the imaginary refractive index 
parameter value, and other model uncertainties have a much 
stronger effect on this parameter. 

Proper knowledge of aircraft attitude is essential for accu- 
rate RSP retrieval. Aircraft attitude was provided by an Ap- 
planix POS-AV, which combines Global Positioning System 
technology with an inertial monitoring unit to determine po- 
sition, velocity, pitch, roll and heading. There are of course 
limitations to the accuracy of such observations, and differ- 
ences between the orientation of the POS-AV and the RSP 
instrument. To account for these errors, we modify our in- 
formation about aircraft geometry using various features in 
the observed data. For example, (Ru) should be minimal 
when defined with respect to the scattering plane for single 
scattering, so slight modifications to the aircraft pitch and 
yaw are made until this is the case. To understand the ef- 
fect of geometric error, we perturbed the aircraft pitch and 
yaw angles by (0.5°). This choice was based on the sensi- 
tivity of our geometry correction techniques, but is proba- 
bly an overestimate. The official POS-AV pitch accuracy is 
(0.03°), while yaw accuracy is (0.1°) (roll uncertainty is not 
accounted for in this work). Table 4 shows that uncertainty 
in aircraft yaw angle is not particularly important (although 
(3 (m f )) uncertainty is large, as it is for many other model 
assumptions). Pitch angle uncertainty, however, has a signif- 
icant impact on the ability to retrieve the fine mode aerosol 
effective variance, cloud droplet effective variance, and es- 
pecially the cloud droplet effective radius. In fact, error in 
the cloud droplet effective radius is larger than the retrieved 
error for that parameter. Accurate retrievals of cloud droplet 
radius thus require accurate monitoring of instrument pitch 
angles. 
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Table 4. This table contains the parameter sensitivity to changes in the forward radiative transfer model. This sensitivity is expressed as 
a percentage of the retrieved parameter value ((100xe/x)). All of the parameter errors from model uncertainties are within the accuracy 
requirements of Mishchenko et al. (2004), and with one exception (the pitch angle uncertainty impact on the cloud effective radius, noted in 
bold), they are all less than the individual retrieval uncertainties. 



t/(0.555) 

3(/w f ) 

r e,f 

v e ,f 

r e,c\ 

u e,cl 

Aerosol layer top raised 30 m 

0.1% 

4.5% 

0.1% 

0.9% 

0.0% 

0.0% 

Cloud layer top raised 30 m 

0.8% 

20.5 % 

0.4% 

4.2% 

0.0% 

0.1% 

Yaw angle changed by (0.5°) 

0.6% 

153.5% 

0.3% 

3.7% 

0.8% 

1.4% 

Pitch angle changed by (0.5°) 

1.0% 

451.5% 

1.6% 

15.4 % 

3.9% 

6.0% 

91 (m f) increased by 0.02 

6.9% 

57.8 % 

2.6% 

12.1 % 

0.0% 

0.2% 

r c (0.555) increased by 0.02 

4.0% 

1481.6% 

0.8% 

16.3 % 

0.0% 

3.3% 

9t(w c ) increased by 0.02 

0.1% 

12.6% 

0.0% 

0.0% 

0.0% 

1.1 % 

k c increased by 0.006 

0.9% 

9.1% 

0.0% 

1.6% 

0.1% 

1.7% 

r e c increased by 0.802 

0.3% 

33.2% 

0.0% 

1.1 % 

0.0% 

0.6% 

Ve^ increased by 0.093 

0.2% 

44.3 % 

0.0% 

1.0% 

0.1% 

0.6% 


As shown in Sect. 3, the information available in RSP ob- 
servations is not great enough to retrieve all optical proper- 
ties for both aerosol modes. We therefore were required to 
assume some of these properties based on prior information. 
In this case, we used the fine mode real refractive index, and 
the coarse mode refractive index and size distribution from 
the Mexico City climatology of Dubovik et al. (2002). We 
used the climatological uncertainty for each assumed param- 
eter to test our retrieval sensitivity. As we can see in Ta- 
ble 4, uncertainty in the assumed fine mode real refractive 
index value has a significant impact on the ability to retrieve 
fine mode optical depth, and a more moderate impact on fine 
mode size distribution. Considering the impact on aerosol 
optical thickness, the retrieval success for that parameter (see 
Fig. 9) indicates that the assumed real refractive index value 
was also valid. Coarse mode optical thickness was deter- 
mined by longest wavelength observations by the AATS sun 
photometer. The assumption is that coarse mode aerosols, 
because of their size, have very little optical thickness spec- 
tral variation. Fine mode aerosols, on the other hand, have 
a very pronounced spectral dependence, and a nearly negli- 
gible optical thickness at long wavelengths. Optical depth 
observations at long wavelengths can therefore be assumed 
to represent the coarse mode only, and that value should be 
constant for all wavelengths. It is difficult to assess the uncer- 
tainty associated with this assumption, but we select a value 
of 0.02, which is larger than the AATS uncertainty alone. Re- 
gardless, coarse mode aerosol uncertainty has a substantial 
effect on fine mode optical depth and effective variance, and 
a large effect on the fine mode imaginary refractive index. It 
appears that the selection of the coarse mode optical depth is 
important, so this must be done with care (or this parameter 
must be retrieved rather than assumed). Uncertainty in as- 
sumptions about the other coarse mode aerosol optical prop- 
erties (complex refractive index and particle size distribution) 
have a minimal effect on retrieved parameter values, and can 
therefore be confidently assumed from climatologies. 


5 Discussion 

This research is a test of the capability of scanning polarime- 
ters to retrieve the optical properties of aerosols lofted above 
clouds (AAC). The test was performed with observations 
from the airborne Research Scanning Polarimeter (RSP). 
The primary strategy for retrieving aerosol and cloud optical 
properties from RSP and APS is to match a multiple scat- 
tering radiative transfer model to observations by iteratively 
modifying model parameters until a match is found. The 
Levenberg-Marquardt method is used to perform this opti- 
mization. The software and tools that we used were devel- 
oped in Knobelspiesse et al. (2011), which couples a pub- 
licly available Levenberg-Marquardt software with a Dou- 
bling and Adding type radiative transfer model. 

An AAC scene observed by the RSP during the Megacity 
Initiative: Local and Global Research Observations (MILA- 
GRO) field campaign was the primary source of data for this 
work. A mixture of urban industrial and agricultural biomass 
burning aerosols from central Mexico were lofted above a 
marine stratocumulus cloud close to shore over the western 
portion of the Gulf of Mexico. Several other instruments, in- 
cluding the Ames Airborne Tracking Sunphotometer (AATS- 
14) and the Solar Spectral Flux Radiometer (SSFR) were on- 
board this aircraft, and collected data that were used to either 
guide the AAC optimization or validate the results. 

This research had three main sections. First, we simulated 
the RSP/ APS retrieval error for an AAC scene with climato- 
logically defined aerosol properties (Dubovik et al., 2002) at 
a variety of optical depths. We found that the aerosol op- 
tical depth, fine mode aerosol size distribution, and cloud 
droplet size distribution can all be retrieved for an AAC 
scene to the accuracy requirements for climate modeling in 
Mishchenko et al. (2004). Aerosol refractive index simulated 
errors are much higher. The error in the real component of 
the refractive index decreases with increasing aerosol optical 
depth, and approach accuracy requirements as optical depth 
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Table 5. Correlation matrix for a simulated AAC scene with (t/(0.555) = 0. 1) and an optimization utilizing angles within (60°) of nadir. 
Elements with absolute values greater than 0.5, indicating strong correlation, are in bold. 



ry (0.555) 

91 (w /) 

S(mf p ) 

%( m f,q) 

r e,f 


t c (0.555) 

r e,cl 

v e,c\ 

z f (0.555) 

1.00 









91 (m f) 

-0.36 

1.00 








S( m.f p ) 

0.26 

0.44 

1.00 







3( mfq ) 

-0.42 

0.27 

-0.55 

1.00 






r e,f4 

0.49 

-0.96 

-0.35 

-0.30 

1.00 





v eJ 

-0.10 

0.64 

0.00 

0.32 

-0.68 

1.00 




r c (0.555) 

-0.40 

-0.45 

-0.53 

0.30 

0.33 

-0.38 

1.00 



A, cl 

0.02 

0.08 

0.07 

-0.01 

-0.01 

-0.01 

-0.00 

1.00 


t’e.cl 

-0.17 

-0.05 

-0.68 

0.30 

0.03 

0.24 

0.01 

0.09 

1.00 


Table 6. Correlation matrix for a simulated AAC scene with (z r (0.555) = 0.1) and an optimization utilizing angles (20°) and (60°) forward 
of nadir. Elements with absolute values greater than 0.5, indicating strong correlation, are in bold. 



t/(0.555) 

91 (m f) 

3 On f p ) 

30 nf.q) 

r e,f 

v e,f 

t c (0.555) 

r e,cl 

v e,c\ 

Tf (0.555) 

1.00 









91 (m f) 

0.44 

1.00 









0.63 

0.18 

1.00 








-0.41 

-0.24 

-0.73 

1.00 






r eJ 

0.10 

-0.82 

0.06 

0.12 

1.00 





v e,f 

0.38 

0.74 

-0.31 

0.16 

-0.51 

1.00 




z c (0.555) 

-0.90 

-0.64 

-0.61 

0.40 

0.15 

-0.42 

1.00 



A, cl 

-0.70 

-0.28 

-0.11 

0.06 

-0.12 

-0.55 

0.49 

1.00 


t’e.cl 

0.71 

0.25 

0.24 

-0.10 

0.19 

0.39 

-0.62 

-0.45 

1.00 


at (0.555 pm) exceeds 0.8. Aerosol single scattering albedo 
error (and the imaginary part of the refractive index from 
which it based) also decrease with increasing optical depth 
to values approaching accuracy requirements. Simulations 
also show that it is preferable to retrieve aerosol properties 
and cloud droplet size distributions simultaneously using all 
available angle observations, rather than to separate the two 
and retrieve cloud droplet properties from observations of the 
cloud bow and aerosol optical properties from observations 
at side scattering angles. We then performed an optimization 
for an example AAC scene. This optimization successfully 
converged to a solution, where a thin layer of aerosols over- 
lay a marine stratocumulus cloud made up of relatively small 
sized droplets. Errors for this optical depth were quite sim- 
ilar to the predictions using simulated data. The retrieved 
aerosols appear similar to the climatology for Mexico City 
aerosols in Dubovik et al. (2002), although with a narrow 
size distribution range possibly associated with particle hu- 
midification or cloud processing. We found evidence of an 
increase in absorption for the shortest wavelengths, which is 
a likely due to Brown Carbon aerosols and was also observed 
by Bergstrom et al. (2010) for the same aerosols. However, 
the imaginary refractive index retrieval uncertainty (which 


determines aerosol absorption) is extremely large. While this 
was predicted with the scene simulations for aerosols with a 
thin optical depth, it limits any conclusions that can be drawn 
from these results. Nonetheless, optical depth and single 
scattering albedo results compare favorably with measure- 
ments by other instruments. Finally, we assessed the impor- 
tance of uncertainty in various model components. We found 
that the aerosol and cloud layer height resolution as provided 
by the CALIPSO instrument will be sufficient for AAC re- 
trievals by scanning polarimeters in orbit. Accurate measure- 
ments of aircraft geometry, especially pitch, are important 
for retrievals with airborne instruments such as the RSP data. 
Finally, while the coarse mode optical depth has a strong ef- 
fect on retrieval uncertainty, assumptions about other coarse 
mode aerosol optical properties only have a small effect on 
the fine mode properties that are being retrieved. 

This work suggests that scanning polarimeters are capable 
of retrieving the total aerosol optical depth, fine mode size 
distribution, and cloud droplet size distribution for aerosols 
lofted above clouds. This is true provided that the scene 
can be simulated as plane layers within the observed pixel, 
and that the layer heights can be supplied by external ob- 
servations (such as the CAFIPSO lidar). Fine mode aerosol 
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refractive index retrievals have a high error that decreases 
with increasing AAC optical depth. These errors approach 
modeling requirements as optical depths approach and ex- 
ceed about 0.8 at a wavelength of (0.555 pm). Most AAC op- 
tical depths will probably not be that large. For lower optical 
depths, we may need to fix refractive index values using a cli- 
matology or external measurements, or use an optimization 
method such as Philips Tikhonov regularization (Hasekamp, 
2010). The retrieved optical depth will be sensitive to the 
accuracy of that assumption. Considering this dependency, 
an alternate approach might be to fix the optical depth (since 
observations of optical depth are much more common than 
those of refractive index), and allow refractive index param- 
eters to vary. We did not attempt this technique in this work, 
but it may be a worthwhile approach in the future, if the RSP 
or APS are observing in coordination with lidars such as the 
High Spectral Resolution Lidar (HSRL), which retrieve both 
aerosol vertical distribution and aerosol optical depth (Hair 
et al., 2001; Hair et al., 2008; Rogers et al., 2009). 

6 Conclusions 

The direct radiative effect of aerosols suspended above 
clouds (AAC) are a potentially important but rarely observed 
phenomenon. Scanning polarimeters, such the Research 
Scanning Polarimeter (RSP), have the capability to retrieve 
AAC properties using nonlinear optimization methods that 
match a radiative transfer model to instrument observations. 
In this paper, we tested the capability of such a retrieval us- 
ing RSP observations of a mixed urban and agricultural fire 
smoke aerosol above a marine stratocumulus cloud in the 
western portion of the Gulf of Mexico. Provided that the ver- 
tical structure of the cloud and aerosol are known, the aerosol 
optical depth, cloud droplet size distribution and aerosol fine 
size mode size distribution can be retrieved. Reasonable val- 
ues were found for the complex refractive index of the fine 
aerosol size mode, but their associated uncertainty is large. 
We also performed a sensitivity study for a similar scene at 
a variety of AAC optical depths, and found that uncertain- 
ties decrease as optical depth increases. If the aerosol opti- 
cal depth exceeds roughly 0.8 at a wavelength of (0.555 pm), 
refractive index uncertainties decrease to a degree that is al- 
most sufficiently accurate for the needs of climate models. 
Once launched, the APS will be able to provide much more 
information about AAC than was previously available. How- 
ever, the complex refractive index will be limited in its accu- 
racy unless external data can be used as a constraint or the 
aerosol optical depth is unusually large. 


Appendix A 

Correlation in simulations 


The analysis in Sect. 3 describes the parameter sensitivity of 
the retrieval, but does not express the correlations between 
the parameters. The off diagonal elements of (C*) (equation 
8) can be used for this purpose if they are normalized by the 
parameter error to create the correlation matrix 


Pij = 


-X,1J 


x/C x Jix/C x ,jj 


(Al) 


Tables 5 and 6 are the correlation matrices for the 
(j f (0.555) = 0.1) simulations in Sect. 3. 
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